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Abstract 

We consider the problem of asymptotic stability of a self-similar attractor for a simple 
semilinear radial wave equation which arises in the study of the Yang-Mills equations in 5-f 1 
dimensions. Our analysis consists of two steps. In the first step we determine the spectrum 
I of linearized perturbations about the attractor using a method of continued fractions. In the 

. second step we demonstrate numerically that the resulting eigensystem provides an accurate 

description of the dynamics of convergence towards the attractor. 

o ■ 

^ ; 1 Introduction 

I ' Self-similar solutions of evolution equations often appear as attractors in a sense that solutions 

\ of an initial value problem starting from generic initial data evolve asymptotically into a self- 

' similar form. In such cases one would like to describe the process of convergence to the self- 

similar solution and understand the mechanism responsible for this phenomenon. These kind 
of problems are relatively well-understood for diffusion equations where the global dissipation 
of energy is the mechanism of convergence to an attractor, however very little is known for 
conservative wave equations where the local dissipation of energy is due to dispersion. In this 



\ paper we report on analytical and numerical studies of this problem for a semilinear radial wave 

equation of the form 

Utt - Urr - -Ur + =y), (1) 

where r is the radial variable, u = u{t,r), and f{u) = — 3u(l — m^). This equation appears 
in the study of the spherically symmetric Yang-Mills equations in 5 -f 1 dimensions (see |21 for 
the derivation). We remark in passing that our results hold for more general nonlinearities, in 
particular for f{u) = sin(2u) which corresponds to the equivariant wave maps from the 3-1-1 
dimensional Minkowski spacetime into the three-sphere. 

It was proved in ^ and later found explicitly in 2 that equation has a self-similar solution 

uit,r)^Uo{p) = ^^^, (2) 

t + gP 

where p ^ r/{T — t) is a similarity variable and T > is a constant (actually, Uq is the ground 
state of a countable family of self-similar solutions [/„ (n = 0, 1, . . . ) but since all n > solutions 
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are unstable they do not appear as attractors for generic initial data). Since 



r=0 



(T - tf 



(3) 



the solution Uq{p) becomes singular at the center when t ^ T . By the finite speed of propagation, 
one can truncate this solution in space to get a smooth solution with compactly supported initial 
data which blows up in finite time. 

In fact, the self-similar solution Uq is not only an explicit example of singularity formation, 
but more importantly it appears as an attractor in the dynamics of generic initial data. We 
conjectured in |3] that generic solutions of equation (Q) starting with sufficiently large initial 
data do blow up in a finite time in the sense that Urr{t, 0) diverges aXt/'T for some T > and 
the asymptotic profile of blowup is universally given by ?7o, that is 



\iu, u{t, iT~t)r)^Uo{r). 

t/ 1 



(4) 



Figure 1 shows the numerical evidence supporting this conjecture. 
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Figure 1: The upper plot shows the late time evolution of some large initial data which blow up at 
time T. As the blowup progresses, the inner solution gradually attains the form of the stable self-similar 
solution Uo{r/{T~t)). The outer solution appears frozen on this timescale. In the lower plot the rescaled 
solutions u{t, (T — t)r) are shown to collapse to the profile Uo{r) (solid line). 
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The goal of this paper is to describe in detail how the limit Q is attained. To this end, in 
section 2 we first study the linear stability of the solution Ua. This leads to an eigenvalue 
problem which is rather unusual from the standpoint of spectral theory of linear operators. We 
solve this problem in section 3 using the method of continued fractions. Then, in section 4 we 
present the numerical evidence that the deviation of the dynamical solution from the self-similar 
attractor is asymptotically well described by the least damped eigenmode. 



2 Linear stability analysis 

The role of the self-similar solution Uq in the evolution depends crucially on its stability with 
respect to small perturbations. In order to analyze this problem it is convenient to define the 
new time coordinate r = — ln(T — t) and rewrite equation in terms of U (r, p) = u{t, r) 

Urr + Ur+2p U pr - (1 - p'W pp + -Up) + ^ = 0. (5) 

In these variables the problem of finite time blowup in converted into the problem of asymptotic 
convergence for t — > oo towards the stationary solution Uq (p) . Following the standard procedure 
we seek solutions of equation in the form U{t,p) = Uq{p) + w{t,p). Neglecting the 0{w^) 
terms we obtain a linear evolution equation for the perturbation w(r, p) 

2 /'(Uq) 

Wrr +Wr +2p Wpr - (1 - p'^){Wpp + -Wp) H -p^ W = 0. (6) 

Substituting w{t, p) ~ e^'^v{p)/p into © we get the eigenvalue equation 

-(1 - p')v" + 2Xpv' + A(A - l)v + = 0, (7) 

where 

We consider equation ^ on the interval < p < 1, which corresponds to the interior of the past 
light cone of the blowup point (t = T,r = 0). Since a solution of the initial value problem for 
equation Q starting from smooth initial data remains smooth for all times t < T, we demand 
the solution v{p) to be analytic at the both endpoints p = (the center) and p = 1 (the past 
light cone). Such a globally analytic solution of the singular boundary value problem can exist 
only for discrete values of the parameter A, hereafter called eigenvalues. In order to find the 
eigenvalues we apply the method of Frobenius. 

The indicial exponents at the regular singular point p = are 3 and —1, hence the solution 
which is analytic at p = has the power series representation 

oo 
n=0 

Since the nearest singularity in the complex p-plane is at p = 1, the series (jSJ is absolutely 
convergent for < p < 1. At the second regular singular point, p — I, the indicial exponents are 
and 1 — A so, as long as A is not an integer, the two linearly independent solutions have the 
power series representations 

oo oc 

vlip) - E ^nHm pT\ v,{p) = bi'\x)ii - p)"+i-\ (10) 

71=0 n=0 
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These series are absolutely convergent for < p < 1. If A is not an integer, only the solution 
Vi{p) is analytic at p = 1. From the theory of linear ordinary differential equations we know 
that the three solutions fo(p)i find V2{p) arc connected on the interval < p < 1 by the 

linear relation^ 

vo{p)^A{\)v^{p) + B{X)v2{p). (11) 

The requirement that the solution which is analytic at p = is also analytic at p = 1 serves as 
the quantization condition for the eigenvalues B(X) — 0. Unfortunately, the explicit expressions 
for the connection coefRcients A{X) and B{X) are not known for equations with more than three 
regular singular points. There are, however, other indirect methods of solving the equation 
B{\) = 0. One of them is a shooting method which goes as follows. One approximates the 
solutions vq{p) and vi{p) by the power series Q and ((Tn|l . respectively, truncated at some 
sufficiently large n, and then computes the Wronskian of these solutions at midpoint, p — 
1/2, say. The zeros of the Wronskian correspond to the eigenvalues (see figure 2). Although 
this method gives the eigenvalues with reasonable accuracy, it is computationally very costly, 
especially for large negative values of A, because the power series © and (tlTill converge very 
slowly. We point out that shooting towards p = 1 fails completely for large negative A because 
the solution V2{p) is subdominant at p = 1, that is, it is negligible with respect to the analytic 
solution vi{p). 
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Figure 2: The Wronskian W[vo,vi]{p = 1/2, A). 



^If 1 — A is a positive integer A'^, then the solution which is analytic at p = 1 behaves as vi ~ (1 — p)"+^ while 
the second solution V2 involves the logarithmic term CjviJi ln(l — p). By a straightforward but tedious calculation 
one can check that the coefficient C]v is nonzero for all A'^. 
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3 The continued fractions method 



In this section we shall solve the eigenvalue problem {Tj) using a method continued fractions. 
The key idea is to determine the analyticity properties of the power series solution vo{p) from 
the asymptotic behavior of the expansion coefficients a„. Substituting the power series Q 
into equation Q we get the four-term recurrence relation (with the initial conditions ao = 1 
(normalization) and a„ = for n < 0) 



where 



P3(n)a„+3 +P2{n)a.n+2 + Pi{n)a„+i + po{n)a„ ^ 0, (12) 



V-i(n) = -100^2 - 850n - 1650, 

P2{n) = -20^2 + (100A-130)n+25A2 + 325A- 750, 

pi (n) = 84^2 + (120A + 378)n + SOA^ + 270A + 618, 

Po W) ^ 36^2 + (36 A + 90)n + 9 A^ + 45 A + 54. 



For n = —2 we have 
and for n — —1 



-350ai + (25A2 + 125A - 570)ao ^ 0, (13) 



-900a2 + (25A^ + 225A - 640)ai + {iOX^ + 150A + 324)ao = 0. (14) 

The series © is absolutely convergent for < p < 1 and in general is divergent for p > 1. In 
order to determine the analyticity properties of the solution Vq{p) at p = 1 we need to find the 
large n behavior of the expansion coefhcients a„. The four-term recurrence relation H12|l can 
be viewed as the third order difference equation so it has three linearly independent asymptotic 
solutions for n — » oo. Following standard methods (see, for example, 0]) we find 

oo (1) / Q\" °° ^(2) / -^X" °° ^(3) 



s=0 \ -/ 



Thus, in general, the solution of the recurrence relation l|12|) behaves asymptotically as 

an - ci(A)a« + C2(A)ai2) + C3(A)a(f). (16) 



If the coefficient ci(A) is nonzero then 



O-n+l "ra+1 
tin 



a„ 



^, (17) 



hence the power series © is divergent for p > 1 (in fact it has a branch point singularity at p = 1). 
On the other hand, if ci(A) = then the solution vo{p) can be continued analytically through 
p = 1. The advantage of replacing the quantization condition -B(A) = in the connection formula 
(jll|l by the equivalent condition ci(A) — follows from the fact that ci(A) is the coefficient of 
the dominant solution in 1)16(1 . in contrast to B{X) which is the coefficient of the subdominant 
solution in 1)11(1 . 

Now, we shall find the zeros of the coefficient ci(A) using the method of continued fractions. 
This method is based on an intimate relationship between three-term recurrence relations and 
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continued fractions. It goes as follows. Suppose that we have a three-term recurrence relation 
(a second order difference equation) 



6n+2 + A„6„+i + B„6„ = 0. (18) 
Let r„ = bn+i/bn. Then, from l|18|l we have 

Br, 



An + r„+i 

and applying this formula repeatedly we get the continued fraction representation of r„ 

Bn Bn+l Bn+2 



An An 



(19) 



A theorem due to Pincherle 4 says that the continued fraction on the right hand side of equation 
l)19|l converges if and only if the recurrence relation H18|l has a minimal solution 6™™, i.e. the 
solution such that lim„^oo &™"V^n — for any other solution &„. Moreover, in the case of 
convergence, equation 119|) holds with r„ = fc™"/^™*" fo^' each n. 

We cannot yet apply this theorem because our recurrence relation H12|) has four terms. How- 
ever, let us observe that 



exact ' 

a„ 
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1 + ^(A-1) 



(20) 



is the exact solution of our four-term recurrence relation (|12|l (although it does not satisfy the 
initial conditions). Using this solution we can reduce the order by the substitution 



3n + 2+ ^(A- 1) 



(21) 



5n + l + ^(A-l) 
to get the three-term recurrence relation 

q2{n)bn+2 + gi(n)6„+i + qo{n)bn = 0, (22) 
where (using the abbreviation 7 = 5(A — 1)/16) 



, , , , 3n-h3-t-7 , , /3\^n-h4-h7 , , 

qo(n) = pi(n) - - — — — — P2(n) + [ - \ — , „ , P3[n), 
5n + 2 + j \5 J n-|-2-|-7 

371 + 4 + 7 

q2{n) = Psin). 

The two linearly independent asymptotic solutions of the recurrence relation (|22|l are 

00 „(1) / q\" °° /3(2) 

a^:^ ^ 4] n-^y/-l^, (23) 



s=0 ^ ^ s=0 

SO in general 

&„^Ci(A)6W+c'2(A)&12). (24) 

Now, our quantization condition ci(A) = is equivalent to Ci(A) = which is nothing else 
but the condition for the existence of a minimal solution for equation H22(l . Thus, we can use 
Pincherle's theorem to find the eigenvalues. 
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In our case An ~ qi{n) / q2{n) and _B„ = qo{n)/q2{n). From ((T^ and we get 
'25A2 + 125A- 570 96 + 15(A - 1) 



bn = 



350 80 + 25(A-l), 

Using Pincherle's theorem and setting n = —1 in l|19|) we obtain the eigenvalue equation 
bo _ 25A2 + 125A- 570 , 96 + 15(A - 1) _ S_i(A) Bo{X) Si (A) 



350 



80 + 25(A-l) A_i(A)- Ao(A)- Ai(A)- 



(25) 



(26) 



The continued fraction in H26|l . which by Pincherle's theorem is convergent for any A, can be 
approximated with essentially arbitrary accuracy by downward recursion starting from a suf- 
ficiently large n — N and some (arbitrary) initial value rjv- The roots of the transcendental 
equation H26(l arc then found nmnerically (see table 1). 



n 





1 


2 


3 


4 


5 


A,i 


1 


-0.588904 


-2.181597 


- 3.570756 


- 5.043294 


-6.486835 


n 


6 


7 


8 


9 


10 


11 


An 


-7.912777 


-9.298265 


-9.907103 


-10.792456 


-12.153033 


-13.164487 



Table 1: The first twelve eigenvalues. 

We recall from [2] that the eigenvalue Aq = 1 is due to the freedom of changing the blowup time 
T. Although this eigenvalue is positive, it should not be interpreted as the physical instability 
of the solution Uo - it is an artifact of introducing the similarity variables and does not show up 
in the dynamics for u{t, r). Note that, strangely enough, all the eigenvalues are real. 



4 Convergence to the attractor 

According to the linear stability analysis presented above the convergence of the solution u(i, r) 
towards the self-similar attractor Uq should be described by the formula 

M(t,r) = C/o(p) + ^Cfee^'=^Wfc(p)/p~ C/o(p) + cie^i^wi(p)/p as r ^ cx), (27) 
fe=i 

where Vk{p)/ p are the eigenmodes corresponding to the eigenvalues A^. and Ck are the expansion 
coefficients. In order to verify the formula H27|l we solved equation numerically for large 
initial data leading to blowup, expressed the solution in the similarity variables, and computed 
the deviation from Uq for t T . Figure 4 shows that for small T — t the deviation from Uq is 
very well described by the least damped eigenmode wi, in agreement with the formula (|27|l . For 
larger T — t the contribution of higher modes has to be included (see figure 5). 



7 



0.1 



LO 



9 0.01 



0.001 




Figure 3: In this figure we determine the coefficient ci in the formula l|27^ as follows. Taking the second 
derivative of ^ at r = we obtain (T - tfurr(t, r = 0) + 16/5 = 2ci (T - . For small T - t we 
plot the left hand side of this equation in the log-log scale and fit the parameters ci and Ai. We get 
Ci — —0.031436 and Ai = —0.58803 (in good agreement with table 1). 
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Figure 4: We plot the deviation of the dynamical solution u{t, r) from the self-similar solution Uo at 
some moment of time close to the blowup time. The solid line shows the least damped eigenmode 
ci(T — t)~^^vi{p)/ p with the coefficient ci obtained in figure 3. 
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Figure 5: The same plot as in figure 4 but at an earlier moment of time. The solid line showing the 
least damped eigenmode (with the same coefficient ci as before) does not fit well the numerical data 
because the other modes have not yet decayed. Including the suitably normalized second eigenfunction 
in the expansion 1)27^ we get a much better fit (dashed line). 
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